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Three bodies moving in a periodic orbit under the influence of Newtonian gravity ought to emit 
gravitational waves. We have calculated the gravitational radiation quadrupolar waveforms and the 
corresponding luminosities for the 13+11 recently discovered three-body periodic orbits in Newto¬ 
nian gravity. These waves clearly allow one to distinguish between their sources: all 13+11 orbits 
have different waveforms and their luminosities (evaluated at the same orbit energy and body mass) 
vary by up to 13 orders of magnitude in the mean, and up to 20 orders of magnitude for the peak 
values. 
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Direct detection of gravitational waves 00 ought to 
come about in the foreseeable future, due to the sub¬ 
stantial effort made at the operational and/or pending 
detectors. One of the most promising candidates for as- 
trophysical sources of gravitational waves are the coalesc¬ 
ing, Lem inspiraling and finally merging binary compact 
stars 0, [I)]. Binary coalescence is the only source for 
which there is a clear prediction of the signal and an 
estimate of the detection distance limit, as general rela¬ 
tivists have completed numerical simulations of mergers 
of compact binaries, such as neutron stars and/or black 
holes, Refs. 0-0. 

Slowly changing, quasiperiodic two-body orbits are 
weak sources of gravitational radiation, Refs. 3 
only accelerated collapse leads to an increase in energy 
loss. The major part of the emitted energy in a binary 
coalescence comes from the final merger of two neutron 
stars, or black holes, that produces an intense burst of 
gravitational radiation. Of course, such mergers are one- 
off events, never to be repeated in the same system, so 
their detection is subject to their (poorly known) distri¬ 
bution in our Galaxy. It is therefore interesting to look 
for periodic sources of intense gravitational radiation. 

There is now a growing interest in three-body systems 
as astrophysical sources of gravitational waves, Refs. US¬ 
UI . These early works did not find a substantial increase 
in the luminosity (emitted power) from representative 
three-body orbits belonging to three families that were 
known at the time, Refs. fl3l - [22l |. over the luminosity 
from a comparable periodic two-body system [2?| . The 
luminosity of a (quadrupolar) gravitational wave is pro¬ 
portional to the square of the third time derivative of the 
quadrupole moment, see Refs. 00 , which, in turn, is 
sensitive to close approaches of two bodies in a periodic 
orbit [28]. Thus, getting as close as possible to a two- 
body collision without actually being involved in one, is 
a desirable property of the radiating system. 

Recently 13 new distinct periodic orbits belonging to 
12 (new) families have been discovered in Ref. 1231. a s 
well as 11 “satellite orbits” in the figure-eight family (24|. 


Some of these three-body orbits pass very close to bi¬ 
nary collisions and yet avoid them, so they are natural 
candidates for periodic sources of intense gravitational 
radiation. 

In this Letter we present our calculations of quadrupo¬ 
lar waveforms, Fig. [I] and of luminosities, see TablejUand 
Fig. [2]of gravitational radiation emitted by the 13+11 re¬ 
cently discovered periodic three-body gravitating orbits, 
Refs. 0, H3|. We have also calculated waveforms of all 
published Broucke-Hadjidemetriou-Henon (BHH) orbits 
14 +20], which we omit from this Letter for the sake of 
brevity, and because they are closely related to Henon’s 
“criss-cross” one, studied in Ref. 0- The waves of 
the 13+11 new orbits show clear distinctions in form and 
luminosity, thus ensuring that they would be distinguish¬ 
able (provided their signals are strong enough to be de¬ 
tected). 

We consider systems of three equal massive particles 
moving periodically in a plane under the influence of 
Newtonian gravity. The quadrupole moment Iij of three 
bodies with equal masses m n = m, (n = 1,2,3) is ex¬ 
pressed as Iij = X^ n =i m x n x n > where x l n is the location 
of nth body, and the spatial dimension indices i and j run 
from 1 to 3 (with x 1 = x, x 2 = y, x 3 = z). The reduced 
quadrupole Qij is defined as Qij = Iij — | Sij Y^t=i -ffcfc- 
The gravitational waveforms denoted by tifT are, asymp¬ 
totically, 


tt _ 2G* d 2 Qjj ( 1 \ 
rc 4 dt 2 \r 2 ) ’ 


(1) 


where r is the distance from the source, Refs. 0,0 ■ Here, 
TT means (i) transverse (Xn=i hj^n 1 = 0) and (ii) trace¬ 
less (Xn=i hu T = 0): where hi denotes the unit vector of 
the gravitational wave’s direction of propagation. The 
two independent waveforms x of a quadrupolar grav¬ 
itational wave propagating along the z axis, Refs. [8;, 0 
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TABLE I: Initial conditions and periods of three-body orbits. £i(0), i/i(0) are the first particle’s initial velocities in the x and 
y directions, respectively, T is the period of the (rescaled) orbit to normalized energy E = —1/2, 0 is the rotation angle (in 
radians) and (P) is the mean luminosity (power) of the waves emitted during one period. Other two particles’ initial conditions 
are specified by these two parameters, as follows: ®i(0) = — 2 : 2 ( 0 ) = —A, 2 : 3 ( 0 ) = 0, j/i (0) = 1 / 2 ( 0 ) = 1 / 3 ( 0 ) = 0, 2 : 2 ( 0 ) = cfci (0), 
£ 3 ( 0 ) = —2£i(0), 1 / 2 ( 0 ) = t/i (0), 1 / 3 ( 0 ) = —2i/i(0). The Newtonian coupling constant G is taken as G = 1 and the masses are 
equal mi , 2,3 = 1 . 


Name 

ii(0) 

yi(0) 

A 


T 

0(rad) 

(P) 

Moore’s figure eight 

0.216 

343 

0.332 029 

2.574 

29 

26.128 

0.245 57 

1.35 

X 

10 u 

Simo’s figure eight 

0.211 

139 

0.333 568 

2.583 

87 

26.127 

0.277 32 

1.36 

X 

10° 

(M8) 7 

0.147 

262 

0.297 709 

3.008 

60 

182.873 

0.269 21 

2.46 

X 

10° 

I.A.l butterfly I 

0.147 

307 

0.060 243 

4.340 

39 

56.378 

0.034 78 

1.35 

X 

10 5 

I.A.2 butterfly II 

0.196 

076 

0.048 69 

4.016 

39 

56.375 

0.066 21 

5.52 

X 

10 6 

I.A.3 bumblebee 

0.111 

581 

0.355 

545 

2.727 

51 

286.192 

-1.090 4 

1.01 

X 

10 5 

I.B.l moth I 

0.279 

332 

0.238 

203 

2.764 

56 

68.464 

0.899 49 

5.25 

X 

10 2 

I.B.2 moth II 

0.271 

747 

0.280 

288 

2.611 

72 

121.006 

1.138 78 

1.87 

X 

10 3 

I.B.3 butterfly III 

0.211 

210 

0.119 

761 

3.693 

54 

98.435 

0.170 35 

3.53 

X 

10 5 

I.B.4 moth III 

0.212 

259 

0.208 

893 

3.263 

41 

152.330 

0.503 01 

7.48 

X 

10 5 

I.B.5 goggles 

0.037 

785 

0.058 

010 

4.860 

23 

112.129 

-0.406 17 

1.33 

X 

10 4 

I.B.6 butterfly IV 

0.170 

296 

0.038 

591 

4.226 

76 

690.632 

0.038 484 

1.23 

X 

10 13 

I.B.7 dragonfly 

0.047 

479 

0.346 

935 

2.880 

67 

104.005 

-0.406 199 

1.25 

X 

10 6 

II.B.l yarn 

0.361 

396 

0.225 

728 

2.393 

07 

205.469 

-1.015 61 

2.33 

X 

10 6 

II.C.2a yin-yang I 

0.304 

003 

0.180 

257 

2.858 

02 

83.727 

0.659 242 

1.31 

X 

10 5 

II.C.2b yin-yang I 

0.143 

554 

0.166 

156 

3.878 

10 

83.727 

-0.020 338 

1.31 

X 

10 5 

II.C.3a yin-yang II 

0.229 

355 

0.181 

764 

3.302 

84 

334.877 

0.472 891 

7.19 

X 

10 10 

II.C.3b yin-yang II 

0.227 

451 

0.170 

639 

3.366 

76 

334.872 

0.254 995 

7.19 

X 

10 10 


can be expressed as 


h+ 

2 g 3 ' 

= -A-'V' rrii(xi 2 +XiXi-yi 2 -yiHi), 

c^r z — ' 

i= 1 

(2) 

h x 

2 Q 3 ' 

= -J-'V' m i (x i y i + 2x i y i + x i yi), 
c^r +—' 

(3) 


i=l 


where r denotes the distance from the source to the ob¬ 
server. We set the units ofG = c = m = l throughout 
this Letter. 

Here the coordinate axes x and y are chosen so that 
they coincide with the orbits’ two (reflection) symmetry 
axes, when they exist, i.e., when the orbits are from class 
I, as defined in Ref. 0. Otherwise, e.g., when only a 
single point reflection symmetry exists, as in class II or¬ 
bits, the x, y axes are taken to be the eigenvectors of the 
moment-of-inertia tensor. The rotation angle necessary 
for each orbit to be aligned with these two axes is given 
in Table HHI. 

The first gravitational radiation waveforms forperiodic 
three-body systems were studied in Refs. |KTll2l . They 
calculated the quadrupole radiation waveforms for three 
periodic orbits of the following three-equal-mass systems: 
(i) of the Lagrange “equilateral triangle” orbit [lj|, (ii) 
of Henon’s “criss-cross” (13: and (id) of Moore’s “figure 
eight” [2lJ. These three orbits are characteristic repre¬ 
sentatives of the (only) three families of periodic three- 
body orbits known at the time. Reference 0 found dis¬ 
tinct gravitational waveforms for each of the three fami¬ 
lies, thus suggesting that one might be able to distinguish 


between different three-body systems as sources of grav¬ 
ity waves by looking at their waveforms [30[ . 

In the meantime 13+11 new orbits belonging to 12 new 
families have been found, Refs. 0 , 0 ], The families of 
three-body orbits can be characterized by their topolog¬ 
ical properties viz. the conjugacy classes of the funda¬ 
mental group, in this case, the free group on two letters 
(a, b), Ref. 0 . The free group element tells us the num¬ 
ber of times the system’s trajectory on the shape sphere 
passes around one or another (prechosen) two-body col¬ 
lision point within one period. Every time the system 
is close to a two-body collision the (relative) velocities, 
accelerations, and the third derivatives of relative coor¬ 
dinates increase, so that the luminosity of gravitational 
radiation also increases; i.e., there is a burst of gravita¬ 
tional radiation. This argument can be made more quan¬ 
titative by appealing to two-body results of Ref. 0 , as 
is shown in footnote [32]. 

We show the gravitational radiation waveforms /i+, x 
in Fig. [D emitted by three massive bodies moving ac¬ 
cording to the orbits from Refs. 0 , 0 belonging to 
these families, where Eqs. © and m are used as the 
definitions of the two waveforms. 


First, we note that all of the calculated three-body 
orbits’ waveforms are distinct [ 31 ]: thus answering (in 
the positive) the question about their distinguishability 
posed in Ref. 0- In Fig. □ we also show the gravita¬ 
tional waveform of one “old” orbit: Simo’s figure eight, 
(discovered in 2002) belonging to the figure-eight fam¬ 
ily. Simo’s figure eight is an important example, as it 
is virtually indistinguishable from Moore’s one, and yet 
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0 T/4 T/2 3T/4 T 


FIG. 1: The gravitational radiation quadrupolar waveforms 
h+ : x x r as functions of the elapsed time t in units of the period 
T, for two periodic three-body orbits (in units of Gm/c 2 ; we 
have set G = m = c = 1 throughout this Letter) and r is the 
radial distance from the source to the observer. Dotted (blue) 
and solid (red) curves denote the + and x modes, respectively. 
Top: Simo’s figure eight, Ref. [22j]; and bottom: orbit I.B.l 
Moth I. Note the symmetry of these two graphs under the 
(time-)reflection about the orbits’ midpoint T/2 during one 
period T. 


the two have distinct gravitational waveforms, see our 
Fig. |T|and Fig. 2 in Ref. |I(ij. That is so because these 
two figure-eight solutions have distinct time dependences 
of the hyperradius R, where R 2 ~ (1 /rn)$ij Ikk, so 
that the two orbits have different quadrupolar waveforms. 

Note, moreover, the symmetry of the waveforms in Fig. 
□ with respect to reflections of time about the midpoint 
of the period T/2: this is a consequence of the special 
subset of initial conditions (vanishing angular momentum 
and passage through the Euler point on the shape sphere) 
that we used. There are periodic three-body orbits, such 
as those from the BHH family, that do not have this 
symmetry. 

The gravitational waveforms’ maxima range from 20 to 
50 000 in our units, with the energy fixed at E = —1/2. 
This large range of maximal amplitudes is due to the 
differences in the proximity of the approach to two-body 
collisions in the corresponding orbits. One can explicitly 
check that the bursts of gravitational radiation during 
one period correspond to close two-body approaches. 

As stated above, the (negative) mean power loss 
( dE/dt ) of the three-body system, or the (posi¬ 
tive) mean luminosity (emitted power) of quadrupo¬ 
lar gravitational radiation (P), averaged over one pe¬ 
riod, is proportional to the square of the third time 

/o\ 

derivative of the (reduced) quadrupole moment Q y - k , 
(dE/dt) = -<P> = -^(G/c 5 )E-, fe=1 (QS ) Q jfe ) = 

~|(G/c 5 ) J2^,k=i (QfkQfk), ( for an original derivation 

see Refs. Si, for pedagogical ones, see Refs. Si). 
(3) 

But, Qj k are proportional to the first time derivatives 


of the gravitational waveforms = {d/dt)Q^ k oc 

(d/dt)/i+ iX . The peak amplitudes of gravitational wave¬ 
forms ft+, x , io turn, grow in the vicinity of two-body 
collisions [321 ]. which explains the burst of gravitational 
radiation as one approaches a two-body collision point. 

The mean and instantaneous luminosities, expressed 
in our units, of these orbits, normalized to E = —1/2, 
are shown in Tableland Fig. (2[ respectively. Note that 
in Table U we show only three of the 11 orbits belonging 
to the figure-eight family: Moore’s, Simo’s, and the sta¬ 
ble choreography (M8) 7 ; they have all the same order of 
magnitude of the mean luminosity [33| , whereas the but¬ 
terfly I and butterfly II orbits, which belong to the same 
topological family, have mean luminosities that differ by 
more than a factor of 40. 

Generally, the mean luminosities of these 24 orbits 
cover 13 orders of magnitude, ranging from 1.35 (Moore’s 
figure eight) to 1.23 x 10 13 (I.B.6 butterfly IV) in our 
units; see Table [T] The peak instantaneous luminosities 
have an even larger range: 20 orders of magnitude; see 
Fig. [2] Here, the symmetric form of the instantaneous 
(time unaveraged) power P = ^(G/c 5 )J2'j,k=iQfkQfk 
was used. This gives us hope that at least some of 
these three-body periodic orbits can, perhaps, lead to 
detectable gravitational radiation signals. 

It is a different question if some or all of these sources 
of gravitational radiation would be observable by the 
present-day and the soon-to-be-built gravitational wave 
detectors: that strongly depends on the absolute values 
of the masses, velocities, and the average distances be¬ 
tween the three celestial bodies involved, as well as on 
the distribution of such sources in our Galaxy. 

Moreover, note that all of the newly found and an¬ 
alyzed three-body orbits have zero angular momentum, 
and many of them are unstable. It is well known H- 
[§0j that by changing the angular-momentum within the 
same family of three-body orbits, the stability of an or¬ 
bit changes as well. So, it may happen that a previously 
stable orbit turns into an unstable one, and vice versa. 
For this reason it should be clear that a careful study 
of gravitational-radiation-induced energy- and angular- 
momentum dissipation is necessary for these orbits [34|. 
Moreover, if realistic results are to be obtained, post- 
Newtonian approximations will have to be applied in the 
future. Such relativistic corrections are most important 
at large velocities, i.e., precisely near close approaches 
that are so crucial for large gravitational radiation. Thus, 
the present Letter is meant only to highlight the possibil¬ 
ities in this field, and should be viewed as an invitation 
to join in the more realistic future studies. 
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FIG. 2: The instantaneous (time unaveraged) luminosity P of quadrupolar gravitational radiation emitted from periodic three- 
body orbits as a function of the elapsed time t in units of the period T. Note the logarithmic scale for the luminosity P (y 
axis). Top left: Moore’s figure eight; second from top left: I.A.2 butterfly II; third from top left: II.B.7 dragonfly; bottom 
left: I.B.l moth I; top right: (M8) 7 ; second from top right: I.A.3 bumblebee; second from bottom right: I.B.5 goggles; bottom 
right: II.B.6 butterfly IV. 


Serbian Ministry for Youth and Sport. 
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